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Abstract 

Over the last decade, the Imaging Air Cerenkov technique has proven itself to 
be an extremely powerful means to study very energetic gamma-radiation from a 
number of astrophysical sources in a regime which is not practically accessible to 
satellite-based instruments. The further development of this approach in recent years 
has generally concentrated on increasing the density of camera pixels, increasing 
the mirror area and using multiple telescopes. Here we present a practical method 
to substantially improve the sensitivity of Atmospheric Cerenkov Telescopes using 
wide-field cameras with a relatively course density of photomultiplier tubes. The 
2-telescope design considered here is predicted to be more than ~3 times more 
sensitive than existing/planned arrays in the regime above 300 GeV for continuously 
emitting sources; up to ~10 times more sensitive for hour-scale emission (relevant 
for episodic sources, such as AGN); significantly more sensitive in the regime above 
10 TeV; and possessing a sky coverage which is roughly an order of magnitude larger 
than existing instruments. It should be possible to extend this approach for even 
further improvement in sensitivity and sky coverage. 
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1 Introduction 



The Imaging Atmosplieric Cerenkov Telescope (ACT) lias proven to be a 
highly successful technique in the detection of high energy, gamma-rays from 
astrophysical sources above ~250 GeV. The catalogue of sources so far de- 
tected by this ground-based technique includes mainly supernova remnants 
(SNR) (see e.g. (U, (0), Q)) and active galaxies (AGN) (see e.g. 0) and ref- 
erences therein), although recently, new classes of objects such as OB star 
associafous (5) aud binary systems fi. have been confirmed to be sources 
of high energy gamma-rays. Also, a recent survey of the galacticplane has 
yielded a significant signal from the region of the galactic center (|7|) and the 
discovery of a series of objects lying along the galactic plane (jH). The nature 
of the physics addressed via high energy gamma-ray measurements includes 
clues to the understanding of the origin of cosmic rays, the nature of AGN 
jets and the density if the extragalactic infrared background. In addition, sev- 
eral fundamental topics of particle physics can be probed, including searches 
for TeV-scale WIMP annihilations (see (0) and references therein), limits on 
radiative neutrino decay |lol ) and a fundamental test of Lorentz invariance 
based on studies of the propagation of very high energy radiation over extra- 
galactic distances (jlll). The sensitivities of such searches are among the very 
best achieved by any technique. 



Almost all current ground-based gamma-ray efforts have concentrated primar- 
ily on improving the sensitivities to lower energies so as to overlap more with 
satellite-based observations and potentially map out a larger range of sources. 
The approaches to this have generally involved increasing the density of cam- 
era pixels, increasing the mirror area and using multiple telescopes. However, 
it is of interest to note that many of the physics topics mentioned above actu- 
ally stand to benefit most from improved measurements at higher energies. It 
is in this regime where AGN spectra and short timescale flux variability most 
critically test acceleration mechanisms as well as uniquely probe the extra- 
galactic infrared background radiation. These higher energy emissions permit 
the direct study of some of the most extreme and poorly understood envi- 
ronments in the universe. There are now several sources known to produce 
gamma-ray emission in the regime above 10 TeV, with evidence of some pro- 
ducing gamma-ray emission up to 40 TeV (see e.g. (0)), and potentially as 
high as 100 TeV (|l^. Establishing the actual range of this emission would, 
in itself, be of substantial interest. Furthermore, it is generally acknowledged 
that much greater sky coverage is desired to study extended sources (j^, to 
better explore the gamma-ray sky and to more continuously monitor multiple 
sources to study transient behaviour. The approaches so far taken have led 
to only a modest improvement in this regard, causing many to suggest the 
need for either a 'brute force' use of many tens of telescopes or some "new 
technology" to allow larger sky coverage at high sensitivity {e.g. ([lih. ([l3)). 
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Our purpose here is to explore a possible method for greatly increasing the 
sensitivity of ACT instruments using an approach which is practical , ec onom- 
ical and does not rely on the development of any new technology (|l5h. This 
approach involves the use of cameras with a relatively large field of view, which 
has the additional advantage of acting as a test-bed for the development of 
"all-sky" instruments. While initially motivated to improve the sensitivity at 
higher energies (as in, e.g., (16)), the approach is seen to lead to improvements 
for energies at least as low as 300 GeV. In this present paper, we will develop 
the basic concept and characteristics of this technique using simulations ap- 
plied to the case of both a single "Wide Field Of View" telescope (WFOV-1) 
and a tandem telescope design (WFOV-2). In a later paper, we will extend 
this approach to a small array and even larger fields of view to explore more 
details of the range and sensitivity of this approach. 



2 The Cerenkov Lateral Distribution of Air Showers and "Wide- 
field" Imaging 

Extensive air showers induced by primary gamma-rays in the regime near 
1 TeV are known to produce a lateral distribution of Cerenkov light at sea 
level with a relatively constant photon density for distances within ~130 m 
from the shower core (0). By a distance of ~200 m from the core, the density 
falls by a factor of 2-3 from this value. Consequently, ACTs have generally con- 
centrated on measuring showers landing within the "inner light pool," where 
the ^constant photon intensity is at its maximum and its relation to the pri- 
mary gamma-ray energy is not particularly dependent on a precise knowledge 
of the core location. This naturally confines the effective areas of individual 
telescopes to the order of a few times lO'* m^. However, this form of the lateral 
distribution is also known to change as a function of both detection altitude 
and primary energy. As an example, figure la shows that, for an altitude of 
2400 m above sea level and for primary energies of ~2 TeV or more, there 
is no longer any distinctive plateau at smaller core distances and the density 
falls in a nearly continuous ~exponential fashion. Therefore, the choice of a 
fiducial radius for analysis becomes somewhat arbitrary. 

If one were to instead base this fiducial radius purely on the observable light 
level, figure la also indicates that the ability of current ACTs to easily trigger 
for energies of less than 1 TeV would imply that a 10 TeV gamma-ray could 
be seen beyond a distance of ~500 m. This would correspond to an effective 
area more than an order of magnitude larger than existing ACTs, and it is 
this attractive possibility which motivates this present work. Owing to simple 
geometry, images from showers at larger core distances tend to have an aver- 
age angle in the camera which is more displaced from the source position. As 
shown in figure lb, an image corresponding to a core distance of 500 m would 
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be displaced by approximately 4 degrees. For point-sources tracked in the cen- 
ter of the field of view, this would require a camera with an angular radius of 
~5 degrees to contain such images. This raises practical concerns regarding 
the need to balance this with a large enough density of photomultiplier tubes 
(PMTs) to adequately define shower image parameters as well as the potential 
impact of mirror aberrations on such images. On the other hand, as indicated 
by figure Ic, shower images become elongated by approximately a factor of 
two at these larger core distances, which should loosen the constraints related 
to these issues to some extent. Another concern is the need to not only main- 
tain the background rejection capability for large image displacements, but to 
actually improve it, since lower intensity images will now have to "compete" 
with a much larger fiux of low energy cosmic-ray background showers (mostly 
due to protons) which bombard the earth from all angles. Finally, there is the 
question of how well the primary gamma-ray energies corresponding to images 
with large displacements can be reconstructed. 



3 Simulated Configuration 

In order to address the issues previously raised, computer simulations were 
used to investigate the "wide-angle" concept in the context of a practical in- 
strument. The simulation package used (^f) was developed based on the EGS4 



checked against various analytic treatments and image characteristics were 
compared both with results from other simulation packages and with actual 
data from the Whipple imaging ACT ('21'). For the purposes of this paper, a 
baseline telescope configuration was chosen which was similar to that of the 
Whipple ACT {e.g. (j22j)), which has a well-documented history of performance 
characteristics for a variety of different camera designs and, thus, allows for 
easy comparison with this device. Accordingly, the model involved a 10 m 
diameter Davies-Cotton mirror design at an altitude of 2400 m above sea level 
and wavelength-dependent mirror reflectivity and PMT quantum efficiencies 
similar to those measured for the Whipple 10 m ACT. However, a slightly 
longer / = 1 focal length was used (compared with / = 0.7), which is more 
in line with current telescope designs and reduces some of the optical aber- 
rations at large image displacements. Each PMT was modelled with a light 
concentrator on the front face, possessing a wavelength-independent reflectiv- 
ity of 80%. The signal recorded in each channel of electronics per event was 
simply assumed to consist of the earliest detected photon arrival time and the 
summed number of photoelectrons (pe) for the corresponding PMT. 



(Hi) and SHOWERSIM 
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Based on the previous discussion, the nominal camera design simulated was 
taken to have an angular radius of 5 degrees, or a full field of view angle of 
10 degrees. This is roughly twice the angular diameter of the cameras used for 
the HESS ACTs fEsh and nearly three times those of the VERITAS project 
(j2J). The cameras for each of these other instruments are composed of an array 
of PMTs which, individually, have a full field of view angle of ~0.15 degrees. 
However, as indicated earlier, gamma-ray images corresponding to shower 
cores which land a few hundred meters away from the ACT are elongated 
by roughly a factor of two relative to those landing within 100 m (see figure 
Ic). Hence, it is reasonable to suspect that PMTs with twice the angular accep- 
tance {i.e. 0.3 degrees) might be appropriate for the approach explored here. 
This would also keep the number of PMTs to less than 1000 per ACT, which 
is certainly in line with current ACT camera design and couches this study in 
very practical terms. Consequently, for the purposes of this paper, the nominal 
simulated camera design consisted of 935 PMTs, each with an angular accep- 
tance of 0.3 degrees, as illustrated in figure 2. Subsequent studies confirmed 
this PMT density to be sufficient and even suggested it could be somewhat 
reduced without substantial impact (section 6). We have chosen to model a 
camera with uniform PMT density (as opposed to one with graded spacing) as 
it both simplifies image analysis and provides a more uniform sensitivity for 
sources across the field of view of the instrument. The minimum event trigger 
conditions required at least any two PMTs in the camera to record a signal 
in excess of 20 pe above fiuctuations in the night sky background, the latter 
of which was scaled from the recorded data of the Whipple 10 m camera and 
corresponded to a Gaussian with an RMS of 3 pe for a PMT with a full field 
of view of 0.25° (15yusr). 



4 Sampling 



In this work, we have concentrated mainly on assessing the performance of the 
hypothetical instrument to point-sources of gamma-rays and have assumed 
such a source to be tracked in the center of the field of view of the cam- 
era. Table 1 summarises the primary simulation data set used in this work. 
For simplicity, gamma-ray events have been generated from the zenith, with 
shower cores randomly sampled over an area on ground within 800 m from 
the detector. Further increase to the effective area for gamma-rays on the or- 
der of a factor of ~3 (though this is energy-dependent) can be achieved when 
viewing sources at larger zenith angles (|25[). but the studies here have been 
limited to what might be achieved relative to the minimum shower 'footprint' 
from overhead sources. The largest core distance used was chosen to ensure 
that the probability of triggering the ACT at the highest gamma-ray energies 
falls to ~zero, so that effective areas can be accurately assessed (figure 3a). 
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Gamma-rays with energies between 300 GeV and 20 TeV were generated based 
on a power-law with a differential spectral index of 2.5, which is indicative of 



energies was also generated for the study of image characteristics at those 
particular energies. Background images were assumed to be predominantly 
due to primary proton interactions in the atmosphere, which constitute the 
majority (~75%) of cosmic-ray showers that trigger ACTs in the TeV regime. 
Heavier nuclei are expected to produce shower images which can be rejected 
at least as easily as those from proton primaries. Images from isolated muons 
are not explicitly considered in this study as these tend to be readily identifi- 
able for the energy regime considered here. The simulated proton events were 
generated from a power-law with a differential spectral index of 2.7 (j28|) and 
sampled over the same area on ground as for gamma-ray showers. To insure 
a reasonable statistical sample in different energy ranges so as to more accu- 
rately characterise the nature of these backgrounds in these regimes, proton 
showers were generated over several, distinct energy intervals with their con- 
tributions appropriately weighted to represent the full spectrum. Cosmic-rays 
are essentially isotropic with respect to the earth, however, in practise, these 
events can only trigger an ACT within a restricted angular range relative to 
the telescope orientation. Accordingly, simulated proton-induced background 
showers were sampled within an angular range of 15 degrees from the telescope 
axis, which was found to be sufficient (see figure 3c). 



To increase the effective statistics, both, gamma-ray and background showers 
were used 100 times each by sampling the Cerenkov photons over different re- 
gions on ground. On average, local variations in the shower properties should 
be relatively uncorrelated. However, within a given shower, all such sampled 
densities are correlated with the first few interactions of the primary particle, 
which determines the height of shower development. Thus, in order to verify 
that the sampling used did not noticeably bias the results of this study, table 
2 compares the average proton and gamma-ray acceptances for different event 
selection criteria (described below) with the case in which repeated sampling 
of each shower was not used. As can be seen, there is no evidence of any 
bias to within the level of the statistical uncertainties of the test. Therefore, 
repeated sampling was used for increased statistics in order to better study 
the differential behaviour of the ACT performance. In cases of low statistics 
in particular parameter regimes, the effect of weighted sampling on the ex- 
tracted differential behaviour was also directly studied to insure there was no 
significant bias due to a handful of heavily-weighted showers. 



many known TeV sources 




smaller sample of 10 fix 
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5 Image Analysis 



For those images wliich triggered the detector, we arbitrarily chose to use the 
same 2-level image cleaning and image parameterisation approach employed 
by the Whipple group for a camera with a similar individual PMT accep- 
tance (0.25 degrees as opposed to 0.3 degrees) (j29|). This required a 4.25(T 
signal above sky noise background for "picture" pixels and 2.25a for neigh- 
bouring "boundary" pixels. In addition, to ensure that enough information 
was available to reconstruct the image parameters reliably, accepted images 
were required to have at least 5 picture pixels. 

Images were then analysed in terms of moments according to the Hillas pre- 
scription (jiO). The relevant image parameters used in this study were: 1) the 
summed number of detected photoelectrons in the image, S; 2) the average 
angular displacement in the camera relative to the source position, V; 3) the 
mean length of the image, as measured along its major axis, C; 4) the mean 
width of the image, as measured along its minor axis, W; and 5) the angle 
of the major axis with respect to the line linking the image centroid with the 
assumed source position, a. All of these parameters were couched in units 
of degrees. Figure 4 shows the distributions of these image parameters for 
gamma-ray and background events. In addition, a further image parameter, 
Ts, was defined as the slope extracted from a linear fit to the earliest pho- 
ton arrival time in each PMT as a function of the angular distance of that 
PMT from the assumed source position. This characterises the projection of 
distances to Cerenkov-emitting regions along the axis of shower development, 
given the expected primary gamma-ray direction. Figure 5 illustrates the def- 
inition of Ts for a particular event. 



6 Dependence on Pixel Size 

The impact of pixel size (PMT density) on ACT sensitivity enters primarily 
though the ability to characterise image shape and orientation so as to aid 
in background rejection. Consequently, we have explored the ability to recon- 
struct the parameters £, W and a as a function of T> for simulated cameras 
corresponding to pixel sizes of 0.22, 0.3, 0.4 and 0.82 degrees in diameter. 
These were also compared with values from actual Whipple and VERITAS 
data for cameras with various pixel sizes over their limited range in V. The 
results are shown in figure 6. Several characteristics are worth noting: 1) The 
ability to determine image orientation (a) is substantially improved at larger 
image displacements. 2) The ability to characterise image parameters is no- 
ticeably improved for pixel sizes of 0.4° compared to 0.8°, but any further 
improvement from using even smaller pixels is relatively modest. This is con- 
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sistent with the findings of Aharonian et al. (j3l|). 3) The fact that very small 
pixel sizes even appear to do worse in determining a at larger image displace- 
ments is an artifact of trigger selection owing to the decreasing photon levels 
in smaller PMTs and the specific trigger requirements assumed. 4) Data from 
actual ACTs seem to have generally larger values of a in the angular dis- 
placement range of ~1° than is predicted by the simulation. Once more, this 
appears to be largely due to the field of view and the fact that smaller cam- 
eras tend to truncate more of the outer parts of these images. This effect is 
illustrated in figure 7, in which comparisons are made with simulations using 
different fields of view. 

In addition, the Point Spread Function (PSF) of the detector due to mirror 
aberrations, has been calculated as a function of the source offset angle (figure 
8). Figure 9 shows the PSF together with the median length {C) and width 
(W) as a function of the angular displacement {V). The figure indicates that 
the PSF should not affect the ability to reconstuction the image angular size. 

Based on the above studies, the choice of the 0.3° pixel size used to explore 
the wide-field concept in this current paper would appear to be more than 
adequate. Indeed, it is likely that a camera with even larger pixels could be 
used without significantly affecting the performance of the instrument. 



7 Single Telescope Performance 

We will first consider a single, isolated, wide-angle telescope and assess its 
performance in terms of background rejection, energy resolution, and effective 
area. Based on this, the sensitivity of the instrument to a hypothetical point- 
source will be estimated for both long term (50 hours) and short term (1 
hour) observation periods, which will be compared with that of other ACT 
instruments. 



7.1 Image Selection Criteria 

In the approach followed here, the goal is to attain good background rejection 
without sacrificing collection area for events over a large range of energies 
and core distances from the ACT. These latter properties are most directly 
related to the parameters S and T). Figures 10 and 11 show the image pa- 
rameter distribution as a function of angular displacement {T>) for gamma-ray 
and background events respectively. Hence, simulated gamma-ray images were 
first categorised into 10 intervals of T) (each spanning 0.5°) and 20 logarithmic 
intervals in »S (each spanning 0.5 in loglO(5)). For each division defined by 
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the resulting 2-D matrix, symmetric acceptance intervals about the median 
of each of the remaining image parameters were defined so as to arbitrarily 
contain 95% of the gamma-ray images. Figure 12 illustrates the definition of 
image selection for events within a given loglO(»S) interval. The intervals for 
the parameters L and W define images which have shapes typical of gamma- 
ray images, whereas the intervals corresponding to a and Ts define images 
whose orientation is characteristic of a primary which came from the assumed 
source position. In particular, a is related to the angular resolution and Ts 
characterises the extent to which the value of T> is due to the displacement 
of the shower core as opposed to a nearby shower with a larger relative incli- 
nation angle which, therefore, is not associated with the source. In addition, 
Ts was required to be positive in order to remove low energy events landing 
very near the detector for which the value of Ts is poorly defined. Table 2 
shows the efficiency of each one of the acceptance cuts individually applied to 
independent simulated gamma-ray and proton spectra. Prom this, it can be 
seen that the introduction of Ts yields an extra factor of between 2 and 3 in 
background rejection. 



7.2 Shower Core and Energy Resolutions 

A minimisation method was employed to simultaneously reconstruct the 
core position and gamma-ray energy of each event so as to explore resolution 
properties of the hypothetical ACT. This approach also has the advantage of 
being trivially extended for use with multiple telescopes, which will be explored 
in section 8. Three image characteristics were used in this fit (see figures 13 
and 14): V, \ogiQ{S) and \ogi(^{Ts/V). The factor of V in the latter expression 
removes the first-order dependence on this parameter and allows the to be 
better approximated by a linear sum of independent contributions from each 
of these characteristics. To define the niodel, simulated gamma-ray events 
were first categorised into 10 m intervals of core distance and logarithmic 
intervals of primary gamma-ray energy using 4 bins per decade. The model 
was extended in energy, with respect to our standard data set (table 1), down 
to 30 GeV and up to 55 TeV to reduce reconstruction biases for events near 
the boundaries of this set. For each division of this 2-D matrix, the median 
value of each image characteristic was determined, along with "equivalent" 
upper and lower one standard deviation values, defined as those values which 
bound ±34% of images about this median. The function was then defined 
as: 



3 r 



rrii — m'^{E, R) 
<Tf{E,R) 



i=l L 
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where rrii is the measured median value of the ith image characteristic for the 
event; and m^^{E,R) and a^{E,R) are the model predictions for the median 
and effective standard deviation for this characteristic, respectively, as a func- 
tion of primary energy, E, and distance from the shower core, R. Appropriate 
cr^{E, R) values were used depending on whether the was above or below 
m'^{E,R) 

Using an independent simulated data set (300 GeV-20 TeV) from that used to 
generate the model (30 GeV-55 TeV), this was then minimised with respect 
to E and R for each event, using values from the 2-D matrix described above. 
Despite the approximations of independent image characteristics and normal 
errors, the resulting distribution of x^, shown in figure 15, is in reasonable 
agreement with the expected 1 degree of freedom. To remove poorly defined 
images, those events with values of the probabilities less than 10% were 
eliminated from the data set as well as those with reconstructed energies 
outside the considered energy range 300 GeV-20 TeV. Together, these cuts 
remove ~27% of the gamma-ray events that passed image cuts and ~66% of 
proton events. Based on this procedure, figure 16 shows the derived core and 
energy resolutions as a function of actual primary gamma-ray energy and true 
core distance from the ACT. For a single telescope, an energy resolution of 
~25% is achieved above 1 TeV, essentially independent of the core distance. 



7.3 Effective Areas 



The "effective area" for a given class of events was computed as the projected 
area on the ground over which simulated core locations are uniformly sampled 
with respect to the ACT, times the fraction of such events which are retained 
in the data sample. Figure 17 shows the derived effective areas for gamma-ray 
and proton showers as a function of the true primary energy for triggered 
showers both before and after image selection. As expected, approximately 
(0.95)^ ~80% of gamma-ray images are retained after selection based on the 
4 image parameters. The effective collection area for gamma-rays after this 
selection is ~ 2 x 10^ m^ at primary energies of 1 TeV and exceeds ~ 5 x 10^ m^ 
above 5 TeV. This is roughly an order of magnitude larger than a conventional, 
single telescope and about twice that of the HESS or VERITAS arrays. At 
higher energies, the relative gain in effective area is even greater. 

For proton-induced showers, the effective collection areas are suppressed by 
roughly a factor of 10^ following image selection. However, it should be noted 

that background rejection factors cannot be directly implied from these figures 
since the true primary energy is not an observable. Thus, the resolution of 
"equivalent" gamma-ray energies and their effect on assumed background and 
source spectra must be taken into consideration. 
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The raw trigger rate from background showers based on an assumed proton 
flux of = 8.73 X IQ-^EtIv m"^ s'^ sr'^ TeV'^ (0), is found to be ~ 

1 kHz, using the basic trigger scheme assumed here. A simple pattern trigger 
can further reduce this by a factor of ~2 (Q, making it easily manageable 
by relatively standard data acquisition electronics (see e.g. (js^). 



7.4 Single Telescope Flux Sensitivity 

The image selection and energy reconstruction procedures previously described 
were applied to simulated proton and gamma-ray spectra sampled from dif- 
ferential power laws of the form E'"^'^ and E~'^'^, respectively. In order to 
couch the derived quantities in terms of typical detection rates, the simu- 
lated spectra were initially normalised according to ^p{E) = 8.73 x lO'^E^^y 
m~^ s^^ sr~^ TeV^^, based on cosmic ray primary measurements ({2^ : and 
J^{E) = 3.2 X IQ-'^E^ev TeV-Vbased on the measured high energy 

flux from the Crab Nebula near 1 TeV (|26|). These were then convolved with 
the relevant effective collection areas and used to yield differential detection 
rates in terms of the reconstructed, effective gamma-ray energy, -E^ec- These 
are shown in figure 18, for different ranges of T>. 

These figures may be used to derive the detector "sensitivity" for such a source 
spectrum by first integrating the rates over a chosen time scale, then scaling the 
assumed gamma-ray flux until a given minimum significance level for detection 
over some defined effective energy range is reached. For the purposes of this 
paper, the sensitivity will be defined based on the integrated signal above 
a given effective energy threshold which yields a detection at a significance 
level of 5(7 and a minimum of 10 detected high-energy gamma-rays. It will 
also be assumed that the nominal background rate is well known. In practise, 
uncertainties in this background rate will reduce the significance level to some 
extent, so that the curves derived here may be more indicative of detection 
levels corresponding to ~3-4o". 

Particularly for this type of camera, background levels will be highly depen- 
dent on both the reconstructed gamma-ray energy and the apparent core dis- 
tance. To account for this, the average signal and background levels were first 
computed separately for the 3 regions of T> shown in figure 18 (0-2.6°, 2.6°-3.7° 
and 3.7°-4.5°) and in each of 8 logarithmic energy intervals between 300 GeV 
and 20 TeV. A maximum likelihood ratio, L^, was then constructed relative to 
the zero-signal case. —2logLj. is approximately distributed as a distribution 
(ji^, from which the average detection significance was calculated assuming 
that three variable parameters of signal strength, power-law index and spec- 
tral cut-off may be fit to the data. This approach appropriately weights the 
contributions from the different signal/background regimes. 
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The resulting sensitivity curves as a function of -E^iec shown in figures 
19 and 20 for exposure times of 50 hours and 1 hour, respectively. These 
are compared with similar curves for the Whipple, MAGIC, VERITAS and 



HESS instruments, taken from (1241 ) and (|35l). To verify the consistency of the 
methods used to obtain the following results, the sensitivity of the Whipple 
telescope was also independently estimated for 50 hours exposure time using 
the same simulation software and similar analysis procedure as employed for 
the wide-angle camera study. As can be seen in figure 19, the results are in 
excellent agreement with published values. 

For continuous sources based on 50 hours of exposure, the wide-angle camera 
is found to be approximately a factor of two more sensitive than the Whipple 
instrument, even at energies as low as 300 GeV. This stems from the fact 
that, even at these low energies, the collection area for the larger field of 
view is roughly a factor of two larger than the Whipple instrument and the 
Ts selection (made more effective by the larger field of view) suppresses the 
background by an addition factor of ~2. As the sensitivity is proportional to 
the square root of both these factors, the overall improvement in sensitivity 
by a factor of two is as expected. 

For this same exposure, the device becomes comparable to VERITAS and 
HESS at energies above 10 TeV, despite backgrounds, owing to the greatly 
increased collection area. For studies based on 1 hour of exposure (as might 
be characteristic of transient emission from AGN), the wide-angle camera 
surpasses the sensitivity of the arrays above 2-3 TeV, despite being only a 
single telescope with a ~30% smaller mirror area than each element of the 
arrays. This is due to the fact that, with it's much larger collection area, the 
wide-angle camera tends to be background dominated rather than statistics 
limited. Thus, any decrease in the relative background level (whether through 
improved background rejection or higher source fluxes over shorter exposure 
times) results in greater improvement relative to conventional cameras, which 
become starved for signal above ~1 TeV (as indicated by the upwards turn in 
their sensitivity curves in figure 20). 



8 Tandem Telescope Configuration 



Background rejection is clearly a key factor to instrumental sensitivity. This is 
particularly true for the wide- field concept, which seeks to extend the dynamic 
range by essentially swapping a statistically limited regime for one which is 
more background limited, as previously mentioned. One obvious way to gain 
further background rejection and also improved energy resolution is with the 
use of multiple telescopes. In order to explore the potential gain in sensitivity 
by this approach, we have limited our consideration here to a tandem config- 
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uration in which identical wide-field ACTs are separated by 125 m, which is 
close to the optimum distance suggested by Aharonian et al. based on trigger 
rates (jsGh. However, note that we have not specifically optimised this configu- 
ration in terms of background rejection versus effective collection area for the 
instrument considered here, but have only chosen this as a "representative" 
configuration to explore the basic characteristics. We have also not optimised 
the trigger for a 2-telescope system, but have simply imposed the same single 
telescope trigger and image analysis/rejection previously described to both 
cameras (figure 21). The simulated data set used for this study is summarised 
in table 3. As with the case of protons in the previous study, showers were 
generated over several different energy intervals and joined via an appropriate 
weighting to represent full spectra so as to improve the statistical characteri- 
sation in very different energy regimes. 



8. 1 Tandem Shower Core and Energy Resolution 

The approach for determining primary energy and core location was mod- 
ified for the 2-telescope case to simultaneously fit both camera images. This 
was done by fitting for the 2-dimensional core location (as opposed to simply 
radius) and introducing a 4th image characteristic, A$, which is sensitive to 
the angular orientation of the core with respect to the coordinates of each tele- 
scope. This parameter was defined as the angular deviation between the core 
direction and the direction of the image centroid, both relative to the centre 
of the particular camera (which is assumed to track the source). This defini- 
tion is illustrated in figure 22. Thus, the function for the "jth" telescope 
becomes: 

rrii — m^(E, R) ^ 
^fiE,R) 

where the sum is over each of the 4 image parameters. The overall function 
for an array of n telescopes is therefore: 

n 

2 _ \ " 2 

^array / ^ 

i=i 

Given the three free fit parameters (the 2-D core location plus primary gamma- 
ray energy), for the case of 2 telescopes (each involving 4 measured charac- 
teristics) , this should nominally correspond to 5 degrees of freedom. However, 
this once more assumes that the image characteristics are completely uncor- 
related and that the uncertainties are normally distributed, neither of which 



=1 
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is strictly true. The actual distribution is compared with the idealised case 
in figure 23. 

In applying this fitting procedure to simulated data, as in the case of a single 
telescope, events with values of T> greater than 4.5° or fit values less than 
10% likely were eliminated from the data set in order to remove poorly defined 
images. Figure 24 shows the derived core and energy resolutions as a function 
of actual primary gamma-ray energy and true core distance from the ACT. 
An energy resolution of ~20% is achieved, as compared to the single telescope 
value of ~25%. 



8.2 Inter- Teles cope Timing and Shower Maximum 

The relative arrival times of the Cerenkov signal in different telescopes can 
also provide valuable information. The thickness of the Cerenkov wavefront 
from an extensive air shower is expected to be ~5-10 ns and the Davies- 
Cotton design considered here results in light paths to the camera which differ 
by as much as 5 ns. However, a typical gamma-ray image involves several 
hundred photoelectrons, such that the statistical sampling should allow a sub- 
nanosecond measurement of the relative wavefront arrival time. For an array 
of telescopes spread over a baseline of ~ 100-200 m, this alone could provide 
a determination of the shower arrival direction with an accuracy comparable 
to or better than that based on image orientation. However, even for the case 
of two telescopes explored here, the relative arrival times of the Cerenkov 
signal can be combined with the reconstructed shower core position to yield 
information on the height of shower maximum and also provide another useful 
discriminant against background showers. 

In particular, if it is assumed that the Cerenkov light is emitted from a height of 
maximum shower development, hmax, above the detector level, for the simple 
case here in which the source is at zenith, the time at which the telescope 
receives the light is then just given by: 

where n is the refractive index of air and R is the distance from that telescope 
to the shower core at the detection level. Thus, in the limit where hmax ^ R 
and n — 1.0, the difference in time between two different telescope measure- 
ments is given by: 

At:^-^{Rl-Rl) 
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Thus, a plot of At versus — Rj) should yield a well-defined hnear relation 
for gamma-rays, whose slope is proportional to the height of shower maximum 
(figure 25). In practise, the lateral extent of the gamma-ray shower at the point 
of maximum development means that much of the light arrives from distances 
somewhat closer to the detector than is implied by the shower core. This 
causes the apparent h„iax derived from the simplistic expression above to be 
overestimated. In order to apply this to a simulated data set, the reconstructed 
core locations were used (as previously described) and the relative "image 
time" in each telescope was defined as that of the image centroid, using the 
linear fit parameters derived in the determination of Ts. Results are shown in 
figure 25 for a simulated gamma-ray source. 

A selection criteria based on this to discriminate against background was 
formed by defining a range of At as a function of {R^ — R\) which contains 
95% of the gamma-rays. This was found to reject ~85% of proton events which 
would otherwise pass data selection. 



8.3 Factorised Background Rejection 

The very large background rejection associated with a multi-telescope system 
can impose statistical limitations on the ability to extract detailed informa- 
tion about the behaviour of the expected background from a finite set of 
simulations. However, many of the image selection criteria are independent 
of each other to a very large extent, thus allowing their individual effects to 
be well determined and then combined as a product. To explicitly verify this 
approach for the current case, we first consider the background rejection of a 
single camera for proton events which have triggered both telescopes and pass 
the selections V < 4.5°, Ts > and npicture > 5. Table 5 shows the fraction 
of such events which pass different data selection criteria, both individually 
and in various combinations. These numbers indicate that selections based on 
1) image shape {£, and W), 2) image orientation a, and 3) hmax, based on 
relative image times (in this case making use of information from the second 
telescope), are all largely independent of each other. 

We now consider these data selection criteria applied to the case of two tele- 
scopes. While the effect of any given selection is generally correlated between 
both cameras, the joint application of each selection will still be independent 
from the joint application of another selection. Table 6 lists the two-telescope 
proton rejection factors for each the three selections listed above. 

Furthermore, the effect of a selection based on the goodness of fit to the 
gamma-ray energy and core location can also be separately applied as a factor 
since the parameters involved (<S, V and Ts) are independent of the three 
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selections above: Image shape and orientation selections are parameterised as 
a function of S and V to explicitly remove any dependence, while the relative 
image times are, by construction, independent of individual telescope time 
measurements. For the case of two telescopes, a selection chosen to keep 
90% of gamma-ray images was found to reject ~85% of proton shower images. 

The combination of these factors for the case of two-telescope analysis is shown 
in table 6 for 2 different ranges of reconstructed gamma-ray energy. As no clear 
dependence on this reconstructed energy is seen, a single, average rejection 
factor will be applied to scale the detection rates from the one telescope to 
the two telescope case. 

8.4 Tandem Telescope Sensitivity 

Following a similar procedure to that employed for a single telescope, figures 
19 and 20 show the sensitivity of a tandem, wide-angle instrument to a hy- 
pothetical point-source with an E'"^'^ differential spectrum, both for 50 hours 
and 1 hour exposure times, respectively. This is compared with similar curves 
for the VERITAS/HESS instruments. The scaling relative to the single tele- 
scope case is consistent with the factor expected from comparing the Whipple 
with the VERITAS/HESS sensitivities Q (figure 19), with an additional 
factor of ~2 gained from the hmax selection. For the 50 hour exposure, the 
sensitivity of the tandem wide-angle design is a factor of ~4 times more sen- 
sitive over the range above 300 GeV, with extended sensitivity above 10 TeV. 
For the case of 1 hour exposure, the sensitivities are comparable at 300 GeV 
and exceed an order of magnitude improvement in sensitivity above 5 TeV. 



9 Summary and Conclusions 

The study presented here suggests that the sensitivity of ACT instruments 
could be substantially improved in the regime above 300 GeV through the use 
of cameras with much larger fields of view for the same PMT density than are 
currently employed. Particularly for the tandem telescope design considered in 
this paper, the ability to image showers further "off-axis" leads to a significant 
gain in effective area with no significant loss in energy resolution compared 
with current instruments, even for showers with cores landing much further 
from the telescopes. The ability to suppress background events also appears 
to be substantial for such a design and yields a very good sensitivity over a 
remarkably large dynamic range. This background rejection has been enhanced 
by the introduction of additional imaging parameters based on both intra- 
and inter- telescope timing characteristics. These parameters may be of some 
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use even for current instruments, though they are particularly well suited for 
cameras with larger fields of view. 

To be explicit, the conclusions of this study rely on four basic premises: 1) 
That the simulations are accurate and that the analysis techniques applied 
are valid. In this regard we have verified the simulation by comparing with 
analytical calculations, alternative simulations and actual experimental data. 
We have been able to reproduce trends seen in other, independent studies and 
replicate parameter distributions and published sensitivity curves from exist- 
ing experiments. 2) That the repeated sampling of showers employed does 
not lead to significant biases. To this end, we have explicitly tested this with 
regard to image selection (Table 2) and have found no evidence of any bias to 
within the limits our statistical uncertainties. Furthermore, throughout this 
analysis, we have specifically checked to insure that resulting distributions 
were not unduly influenced by a handful of independent showers with unusu- 
ally high "weights." 3) That the background rejection factor for the tandem 
WFOV design can be factorised into groups of largely independent contribu- 
tions which can be separately assessed and combined as a product. However, 
in addition for there being logical arguments for why this ought to be the case 
for the parameters used, we have also explicitly verified the independence of 
these parameters to within the limits of our statistical uncertainties (Table 5). 
4) That the predicted rejection factors due to newly introduced timing pa- 
rameters are accurate, even though this is yet to be experimentally tested. We 
find no reason to doubt these factors given that the photon timing is largely 
governed by basic air-shower development and geometry, though we certainly 
encourage future experimental efforts to explicitly explore this. Consequently, 
we believe the basic results of this study to be robust. 

These results predict the 2-telescope design considered here to be more than 
~3 times more sensitive than existing/planned arrays in the regime above 
300 GeV for continuously emitting sources; up to ~10 times more sensitive for 
hour-scale emission (relevant for episodic sources, such as AGN); significantly 
more sensitive in the regime above 10 TeV; and possessing a sky coverage which 
is roughly an order of magnitude larger than existing instruments. This is 
despite the fact that the overall mirror area for the instrument modelled here is 
nearly three times smaller than that of either the current HESS or VERITAS- 
4 designs. This study also suggests that the approach could be extended even 
further by using larger fields of view, mirror areas more comparable to the 
current devices and a more optimised array configuration to both increase the 
effective area further and potentially take better advantage of inter-telescope 
timing so as to improve the angular resolution, even at lower energies. These 
latter issues will be explored in a future paper. 

This research has been supported by the Particle Physics and Astronomy 
Research Council (PPARC). 
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Table 1 



Summary of principle shower simulations used in this study for a single telescope 
(WFOV-1). The numbers in parenthesis are the effective number of events used 
after normalising according to the energy spectrum. 



Selection 

Parau/uM cr 


txi'y 

100 Thix)ws 


1 Throw" 


100 Tlirow-s 1 Throw 


C 


0.852±0.002 


0.84±0.02 


0.174±0.003 


0.18±0.03 


w 


0.852±0.002 


0.85±0.02 


0.222±0.003 


0.28±0.04 


a 


0.875±0.002 


0.87±0.02 


0.037±0.001 


0.03±0.01 


rs 
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0.84±0.02 
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{rs>o 
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0.367±0.005 


0.37±0.04) 


All Cuts 

All but Ts 


0.750±0.002 
0.795±0.002 


0.74±0.02 
0.78±0.02 


0.0017±0.0003 
0.0051±0.0005 


0.0041±0.0041 
0.0044±0.0041 



Table 2 



Effects of image selection as defined in the text, comparing the case where each 
shower is sampled once with that in which each is sampled 100 times, n^^bkg are the 
acceptance values, with respect triggered events, for gamma and background events, 
respectively. Note that only events with 5 or more picture pixels are parameterised 
(which, by itself, removes about 10% of triggered gamma events and about 40% of 
triggered proton events). 
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Table 3 



Summary of principle shower simulations used in this study for the tandem config- 
uration (WFOV-2). The numbers in parenthesis are the effective number of events 
used after normalising according to the energy spectrum. 



Background rejection factors for a single wide-angle telescope. 



Energy (TeV) 


0.3-20 


Image Cut 


Values Relative to Trigger -|- V 


AW 


(7.1±0.4)-10-2 


a 


(1.01±0.05)-10-1 




(l.l±0.1)-10-2 


Product of 




(AW, a) 


(5.3±0.3)10-3 


and 





Tabic 4 



Background rejection factors for a single wide-angle telescope (WFOV-1) and dif- 
ferent combinations of image selection criteria. The cut in distance, V, includes the 
conditions 7^ > and >5 picture pixels in the image. Events are only considered 
as having passed image selection if the Ej^^ is within the considered energy range. 
The cut alone removes ~25% of the background events (see text for details). 
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Background rejection factors for a single 
wide-angle telescope given a two telescope trigger. 



Energy (TeV) 


0.3-20 


Image Cut 


Values Relative to Trigger -|- T> 


AW 


(1.4±0.3)-10-1 


a 


(1.4±0.3)-10-1 


^max 


(2.1±0.4)-10-1 


C,W,a 


(4.9±1.7)-10-2 


-C, W, hmax 


(4.4±1.6)-10-2 




(4.3±1.5)-10-2 


C,W,a, \lmax 


(2.4±1.2)10-2 


Product of 




(£,>V,a, h.max) 


(4.1±1.5)10-3 


and 





Table 5 

Background rejection factors for a single wide-angle telescope and different com- 
binations of image selection criteria. Note that a second telescope has to be used 
to evaluate the trigger and the hmax cut, but that only one telescope is used to 
evaluate the rest of the cuts. The cut in distance, V, includes the conditions Tg > Q 
and >5 picture pixels in the image. Events are only considered as having passed 
image selection if the -E'^gg is within the considered energy range. 
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Background rejection factors for the tandem telescope configuration. 





T2/T1 


Energy (TeV) 


0.3-3 3-20 


0.3-20 
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Values Relative to Trigger -|- V 
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(7.1±0.7)-10-3 
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(3.5±0.2)-10-l 
(7.0±0.8)-10-2 
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Product of 
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(2.8±0.4)-10-6 (4.6±0.8)-10-6 


(3.7±0.4)-10-6 


(6.9±0.4) 10-4 


Relative 

Trigger + D 






(2.9±0.2)10-1 



Table 6 



Background rejection factors for the tandem telescope configuration and different 
combinations of image selection criteria. The cut in distance, V, includes the con- 
ditions Ts > and >5 picture pixels in the image. Events are only considered as 
having passed image selection if the Ej^^,^ is within the considered energy range 
shown in column 1. The cut alone removes ~85% of the background events (see 
text for details). The relative rejection factor of WFOV-2 to WFOV-1 (table 4) is 
shown in the last column of the table (T2/T1). 
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Fig. 1. a) Simulated average Ccrcnkov lateral distribution for showers initiated 
by gamma-rays of various energies, convolved with atmospheric extinction, mirror 
reflectivity and typical photomultiplier tube efficiency. All simulated showers were 
generated from a zenith angle of 0° with an assumed observation level for detection 
of 2400 m above sea level, h) Image angular displacement (P) as a function of core 
distance and c) image length {€) as a function of angular displacement for 0.5, 1, 
2, 10, 20 and 50 TeV gamma-ray showers. The shaded bands include 68% of the 
events about the median. The observed turn-over of the image length distribution for 
higher energies and larger displacements is due to image truncation by the camera 
edge (with a radius of ~5°). 
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Fig. 2. Representative images of simulated 20 TeV gamma-ray showers with in- 
creasing distance parameter, going from a), 0.61° through b), 2.44° c), 3.42° and 
d), 4.41°. Different shaded pixels indicate light intensity. 




a) 



b) 



c) 



Fig. 3. Trigger probability as a function of energy and core distance for simulated 
gamma-ray events a) and background events b), where the solid line and dashed 
line indicate the 10% and 50% trigger probability contours, respectively, in the case 
of gamma events, and 10% and 20% in the case of background events, c) Trigger 
probability for proton showers as a function of zenith angle. 
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Fig. 4. Image parameter distributions for simulated triggering gamma-ray events of 
6 different energies (not normalised according to spectral shape) and background 
events summed over all energies (see table 1). 
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Fig. 5. a) Arrival photon time on the camera for a simulated gamma-ray event. 
Each point in the graph corresponds to the arrival time of the first photon at a 
particular PMT. Only image pixel are used to derive the slope, Ts- The corre- 
sponding gamma-ray image is shown in the in-set figure, h) Distribution of Ts for 
simulated triggering gamma-ray images of different 6 energies (not normalised to 
account for spectral shape). The continuous line is the equivalent distribution of Tg 
for simulated background events summed over all energies (see table 1). 
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Fig. 6. a) Median alpha and b) median length as a function of the image angular 
displacement for cameras with a ~10.° FOV but using 4 different pixel sizes. For 
comparison, data from four different Cerenkov telescope cameras, with different 
pixel sizes and FOVsJiave been added to the plots, Whipple-490 (2I), VERITAS-Tl 
Q, Whipple-109 (Isi) and Whipple-37 Effect of a hmited FOV on alpha 
and d) length. For simplicity, only 1 TeV gamma-ray images are considered on these 
two bottom plots. 
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Fig. 7. Simulated 1 TeV gamma-ray event at dcore=120 m (top panel), and 20 TeV 
gamma-ray event at dcore=120 m (bottom panel) as recorded by cameras with 
different field-of-view (FOV). Prom left to right, ~ 3°, ~ 5° and ~ 10° FOV. The 
first two arc similar in FOV to Whipple cameras and the third one, corresponds to 
the wide field of view camera used in the simulations presented in this work. Note 
as well that the pixel size is different for each camera (from left to right, ~ 0.19°, 
~ 0.23° and ^ 0.3° FOV). Superimposed on top of the shower images are 5% (outer 
line), 10%, 50% and 90% photon density contours. 
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Fig. 8. Distribution of photons on the camera plane for an on-axis point source 
and for two different off-set positions, 1.5° and 4.°. The circle corresponds to the 
individual PMT area. The solid lines in the three bottom plots indicate the expected 
peak position, the dashed lines are the 1 a region along the radial direction ((7^) 
and the dotted lines indicate the PMT size. 
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Fig. 9. a) RMS of the PSF (ur) along the radial direction for different source offset 
positions. The shaded region includes 68% of the events about the median length, C, 
as a function of radial distance for a subset of simulated 2 TeV gamma-ray showers. 
b) Similar plot for the RMS of the PSF along the angular direction (cr^), where the 
shaded region includes 68% of the events about the median width, W, for the same 
showers. 
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Fig. 10. Image parameters as a function of angular displacement on the camera for 
triggering gamma-ray showers of 6 different energies. The shaded regions include 
68% of the events about the median. 
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Fig. 11. Image parameters as a function of angular displacement on the camera 
for triggering proton showers summed over all energies. The shaded regions include 
% of the events about the median. 
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Fig. 12. Acceptance intervals for each of four image parameters used in this work 
for background rejection as a function of angular displacement, P, for a given range 
in Iog^o('5)- The dots correspond to individual gamma-ray images. The width of 
each box indicates the angular displacement interval over which median values are 
obtained (see text for details), while the length of the box corresponds to the 95% 
acceptance interval about the median. 
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Fig. 13. Image size, S, as a function of a) core distance b) angular displacement V 
and c) Ts foT gamma-ray showers of 6 different energies. The shaded regions include 
68% of the events about the median. 
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Fig. 14. Distribution of image centroids on the camera plane for different gamma-ray 
energy showers. 
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Fig. 15. distribution of reconstructed single wide-angle telescope events. Only 
events which energy falls in the 300 GeV-20 TeV energy range are represented here. 
The dashed line is the probability density function for the distribution with 1 
degree of freedom. 
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Fig. 16. RMS energy and Core resolution as a function of true energy (E^^) and true 
core distance (Corer) for a single wide-angle telescope. 
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Fig. 17. a) Gamma-ray collection areas, both at the trigger level and after image 
cuts, for a single wide-angle camera (WFOV-1). The different lines correspond to 
effective areas after image cuts for W hipp le-490 and Whipple-490 at Large Zenith 
Angle (LZA) Hi) (dashed line), HESS Q (continuous hne) and VERITAS-Tl Q 
(dotted line), VERITAS-7 (24 j (dotted-dashed line), h) Gamma-ray acceptance as 
a function of energy, c) Background collection areas, both at the trigger level and 
after image cuts, for a single wide-angle camera (WFOV-1). The collection area 
for protons has been normalised to the solid angle subtended by the camera, d) 
Background acceptance function of energy. 
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Fig. 18. Differential rates for a) gamma-ray and h) background events as a function 
of reconstructed gamma-ray energy, E]^^^. di, d2 and d3 refer to the three angular 
displacement intervals defined in the text. 
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Fig. 19. Integral sensitivity curve for a 50 h exposure for the detectors simulated in 
this work assuming an integral E~^'^ source spectrum. WFOV-1 refers to a single 
detector while WFOV-2 refers to the tandem telescope configuration described in the 
text. Also shown is the sensitivity curve for a Whipple-like detector derived using the 
same simulation software used to explore the wide-angle telescope considered in this 
work. These are compared with the stated sensitivity curves for other experiments 
from (0) and 0. 
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Fig. 20. Scaled 1 h integral sensitivity curves for the detector simulated in this work 
assuming an integral E~-^-^ source spectrum. WFOV-1 refers to a single detector 
while WFOV-2 refers to the tandem telescope configuration described in the text. 
Also shown is the expected response of VERITAS-7 for a 1 h exposure. 
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Fig. 21. a) Comparison of the gamma-ray collection areas (after image cuts) of a sin- 
gle wide-angle camera detector (WFOV-1) and the tandem telescope configuration 
(WFOV-2) simulated on this work, b) Ratio of these effective areas. 




Fig. 22. Illustration of the model parameter A$ used for core reconstruction with 
the tandem telescope configuration. 
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Fig. 23. distribution of reconstructed tandem-telescope events. Only events which 
energy falls in the 300 GcV-20 TcV energy range are represented here. The dashed 
line is the probability density function for the distribution with 5 degrees of 
freedom. 
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Fig. 24. RMS energy and Core resolution as a function of true energy (Ey) and true 
core distance (Corer) for the tandem telescope configuration. 
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Fig. 25. Relative shower arrival time as a function of relative squared distance to the 
shower core to the different telescopes. The dark-shaded band corresponds to the 
95% gamma-ray events which have passed image selection. The light-shaded band 
is the equivalent for background events, where only events which have passed the 
trigger condition and have at least 5 picture pixels and Ts >0. are represented. 
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